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ABSTRACT 


The  theory  of  birefringence  in  a  transparent, 
amorphous  solid  is  developed  from  fundamental 
principles  in  order  to  obtain  a  design  formulation 
for  spatial  light  modulators  that  will  operate  at 
100  MHz  center  frequency  with  at  least  50  P«r  cent 
bandwidth.  Both  longitudinal  and  shtar  wave 
elasto-optical  interactions  are  treated. 

Relations  for  the  phase  modulation  index  and 
the  diffracted  light  intensity  and  polarization, 
as  functions  of  the  modulator  parameters,  are  pre¬ 
sented  for  the  case  of  normal  light  incidence. 
The  effect  of  internal  refraction  on  the  first-or¬ 
der  peak  intensity  is  then  quantitatively  examined. 
Prom  this  consideration,  a  criterion  for  obtaining 
maximum  optical  performance  in  the  modulator  is 
derived. 

An  analysis  of  the  quartz  transducers  used  to 
generate  the  elastic  waves  in  the  modulators  is 
presented  as  an  integral  part  of  the  design  formu¬ 
lation.  The  effect  of  lead  and  indium  bonds  on 
transducer  response  is  treated  in  detail. 
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I.  INTRODUCTION 


Ad  part  of  the  continuing  effort  in  the  development  of 
electro-optical  signal  processing  techniques  for  array  an¬ 
tennas  at  the  Columbia  University  Electronics  Research  Lab¬ 
oratories,  a  research  program  on  solid  light  modulators  hau 
been  initiated  to  extend  the  system  bandwidths  and  time- 
bandwidth  products  of  real-time  processors  beyond  that  which 
is  presently  available. 1-12 

The  program  has  proceeded  along  two  complementary  paths: 
a  detailed  theoretical  analysis  of  birefringence  in  amorphous 
solids,  with  attention  focused  on  its  application  to  spatial 
light  modulators,  and  experimental  verification  of  the  theory 
and  determination  of  the  necessary  empirical  data  for  optimal 
system  design,  so  that,  ultimately,  the  feasibility  of  em¬ 
ploying  these  devices  in  very  high  frequency,  w.ide-band  pro¬ 
cessors  can  be  demonstrated. 

This  report  covers  the  theoretical  work  performed  over 
the  past  year  and  is  concerned,  specifically,  with  the  modu¬ 
lator  design  rather  than  system  considerations.  In  format, 
it  is  divided  essentially  into  two  parts.  In  the  first  part, 
the  theory  of  birefringence  in  amorphous  solids  is  derived 
from  crystal  optics  theory13  and  the  elasto-optic  interaction 
theory  as  developed  by  Mueller.14  In  the  second  part,  the 

t 

fundamental  theory  is  then  applied  to  obtain  a  design  formu¬ 
lation  for  solid  light  modulators  that  will  operate  in  a  CW 
mode  at  80  and  100  MHz  center  frequency. 


1  For  numbered  references  see  Sec,  V. 
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Consideration  has  been  given  to  both  longitudinal  and 
shear  wave  interactions  to  indicate  the  differences  in  modu¬ 
lator  performance  that  can  be  expected,  and  to  establish  the 
electrical  drive  requirements  in  each  case. 

The  treatment  is  restricted  to  fused-silica  solid  light 
modulators  which  are  driven  by  quartz  transducers  that  are 
bonded  to  the  delay  medium  with  thin  bonds  of  lead  or  indium 


s 
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II  theory  of  birefringence  in  amorphous  solids 


When  a  beam  of  light  is  passed  through  a  stressed,  trans¬ 
plant  amorphous  solicit  becomes  birefringent,  i.e.,  it  is 
doubly  refracted.  The  output  light  is  polarized  in  two  spe¬ 
cific  directions  and  the  polarized  components  are  shifted  in 
phase  with  respect  to  the  phase  of  the  incident  light.  The 
amount  of  phase  shift  exhibited  by  the  components,  at  each 
point  in  the  solid,  is  a  direct  function  of  the  stresses  at 
that  point.  Hence,  if  the  stresses  are  zero,  the  phase  shifts 
will  be  zero  and  the  output  light  will  remain  essentially  the 
same  as  the  incident  light.  Thus,  an  amorphous  solid  is  only 
temporarily  double  refracting  to  the  extent  of  the  stress  dis¬ 
tribution  within  it  at  any  instant  of  time. 

Transparent  crystals  such  as  calcite  and  quartz  are  also 
double  refractive.  However,  this  is  an  intrinsic  property  of 
these  solids  and,  therefore,  exists  even  though  they  are  un¬ 
stressed.  The  explanation  for  this  property  is  based  upon  a 
fundamental  assumption  of  the  existence  of  a  permanent  state 
of  electrical  anisotropy  within  the  crystals.  Since  the  char¬ 
acter  of  the  double  refraction  in  crystals  is  similar,  elec- 
tromagnetically ,  to  that  in  stressed  amorphous  solids,  the 
theory  developed  for  the  optics  of  crystals  should  be  appli- 
cabxe  to  the  latter  phenomena,  if  it  is  assumed  that  by  stress- 
ing,  the  normally  electrically  isotropic  solid  is  converted 
into  one  which  is  temporarily  electrically  anisotropic. 

In  the  following  treatment  the  approach  will  be  to  de¬ 
velop  the  electromagnetic  theory  first,  so  that  the  optical 
effects  of  the  alasto-optic  interaction  become  clear.  Once 
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thi.  i.  accomplished,  attantion  will  he  directed  to  elastic 
consideration,  that  cre.t.  an  anisotropic  state  n  ^amor 
P--  solid.  These  elastic  -scions  «  “  ^  u.to_ 
lated  to  the  optic,  in  order  to  derive  the  ru 

optical  equation*. 

-r . ™  A  jsmsmm. 

consider  a  plane,  ^chromatic  light  wave  of 

.  tort  propagating  in  a  transparent  crystal 

”  ttat  is  homogeneous ,  non-conductive,  -9*^ 
isotropic  and . anisotropic.  The  relationship  be¬ 

tween  the  vectorf  amplitudes  of  the  electric  vectors, 

1,  is  then  given  by. 


A. 


Ji 


(k  -  1,2,3) 


(1) 


.  D.  e  are  the  amplitude  components  of  the  field  *ec- 

where  D*,  *  arbitrary  reference  frame  and 

tors  along  the  a.  axes  w,. 

C  are  the  dielectric  constants  of  the  crys 
M  Maxwell 1  s  equations,  under  the  conditions  described  above, 

become: 


V  X  3C  »  ~  0 


(2) 


V  x£ 


-  ±8 
c 


With  the  dot  over  the  vector,  implying  differentiation  with 
respect  to  time,  t. 

If  the  representation  of  the  electromagnetic  field  vec¬ 
tors,  X,e,&,  and  a.  is  of  the  form. 
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*-tl, 

Eq.  (2)  reduces  to 

A  _ 
n*  X  H  >  -  D 

(3) 

ni  X  E  ■  I  ■  yH, 

where  n  is  the  refractive  index  in  the  crystal,  defined  by 


v  is  the  phase  velocity  of  propagation  of  the  light  wave, 

P  A 

c  is  the  speed  of  light  in  vacuum,  s  the  unit  vector  in 

the  direction  of  light  propagation  and  p.  the  magnetic 

permeability  which  will  be  taken  as  equal  to  one. 

Solving  Eq.  (3)  for  D,  in  terms  of  E,  gives  the  fun¬ 
damental  relationship  between  these  vectors  in  the  electri¬ 
cally  anisotropic  solid, 

D  -  n2[E  -  (E  •  1)1).  i'  (5) 

Note  that,  in  general,  D  will  not  be  parallel  to  E  but 
rotated  as  shown  in  Fig.  1.  The  fact  that  D,  E  and  s 

lie  in  the  same  plane  is  a  consequence  of  Eq.  (3),  and  U.1.3. 

» 

The  anisotropy,  as  expressed  by  Eq.  (1),  has  a  geo¬ 
metrical  interpretation  which  is  very  helpful  in  understand- 

* 

ing  and  explaining  the  optical  action  in  crystalline  media. 

To  arrive  at  this  interpretation,  consider  the  expression 
for  the  electric  energy  density. 
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FIG.  1  GEOMETRY  OF  ELECTROMAGNETIC  FIELD  VECTORS 
IN  ANISOTROPIC  MEDIUM 
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3  3 

we  "  Btt  k^1  \ektEf 

where  it  is  assumed  that  wfi  ■  (jjjf)  E  •  D  is  applicable  to 
the  anisotropic  case.  Since  w  is  never  negative,  Eq.  (6) 
is  a  positive  definite  quadratic  form  whose  coefficients 
may  be  shown  to  be  the  components  of  a  symmetric  tensor.  When 
this  condition  exists,  three  mutually  orthogonal  axes  can  al¬ 
ways  be  found  along  which  the  c^Ckj^t)  are  identically  zero. 
Then  Eq.  (6)  reduces  to, 

we  “  far  €k*k2'  ^ 

where  the  are  the  principal  dielectric  constants  of  the 

medium  and  the  E^  are  the  components  of  the  electric  inten- 
sity  in  the  directions  of  the  principal  axes  of  the  quadric 
surface  defined  by  Eq.  (6).  Along  these  principal  axes,  it 
is  seen  that 

D*  -  e^,  (k  -  1,2,3)  (3) 

i.e. ,  is  parallel  to  E^. 

Choosing  the  coordinate  reference  frame  coincident  with 
the  principal  axes  of  the  quadric,  it  is  possible  to  deter¬ 
mine  the  geometrical  character  of  the  surface  by  first  writing 
Eq.  (7)  in  the  form, 

3  °v2 

8ttw  -  k  «  i  ~~  ,  (9) 

e  k-1  *k 

using  Eq.  (8) ,  and  then  defining  new  quantities. 
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1,2,3) 


bo  that  Eq.  (9)  becomes  r 


3  *k 

2  IT 
k-l  S 


2 


1. 


(10) 


(11) 


This  is  recognised  as  the  equation  of  an  ellipsoid  whose 
semi-axes  have  lengths  equal  to  Thus,  the  electrical 

«.nisotropy  of  the  crystalline  medium  at  every  point  is  ex¬ 
pressible  in  terms  of  an  index  ellipsoid  whose  size  and  orien¬ 
tation  is  governed  by  the  dielectric  properties  of  the  medium. 
The  relationship  between  the  index  ellipsoid  and  the  optical 
effects  will  now  be  considered. 


Solving  Eq.  (5)  for  exclusively,  by  applying  Eq.  (8), 

gives 


E*  -  nZ*E  '  *}  «k.  (k  -  1,2.3)  (12) 

n  '  ck 

where  s^  are  the  components  of  s  along  the  principal 
axes  of  the  index  ellipsoid,  taken  as  a  reference  frame. 
Multiplying  each  E^  of  Eq.  (12)  by  its  corresponding  s^ 
and  summing  then  yields, 

E  •  s  «  n2  I  — - -  (E  •  s)  . 

k»l  n2  - 
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Hence , 


3 

Z 

k-1 


1 


n* 


If  this  equation  is  first  multiplied  through  by  n2,  and 
s  •  s  •  1  is  subtracted  from  both  sides,  there  results 


3  n2s.  2 

Z  - — 

k*=l  n2  - 


3  *k 

Z  — - — 
k-1  n2  - 


3 

Z 

k-1 


n 


-  0. 


Applying  Eq.  (4)  to  eliminate  n  in  the  last  expression, 
and  defining  the  principal  velocities  of  propagation  as 


(* 


1,2,3) 


reduces  that  expression  to, 


3 

Z  - - * - 

k-1  v  2  -  v  2 
P  * 


0  , 


(13) 


(14) 


which  is  Fresnel's  equation  of  wave  normals.  This  equation 
permits  the  determination  of  the  phase  velocity  of  the  light 
wave  in  the  medium,  once  the  principal  dielectric  constants 
and  the  direction  of  propagation  of  the  light  are  known. 

In  general,  there  will  be  two  independent  solutions  for 
v  given  by  the  Fresnel  equation.  Hence,  a  linearly  polar¬ 
ized  plane  wave  incident  on  an  electrically  anisotropic  med¬ 
ium  will  be  split  into  two  linearly  polarized  components  that 
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travel  with  different  velocities,  v  and  v_  ,  and  are 

Pi  P  2 

oriented  according  to  the  resultant  of  the  electric  inten¬ 
sity  components  given  by  Eq.  (12),  when  it  is  written  in  the 
fora 


v  2 

E,.  -  — — * - -  (E  •  »  (*  "  1,2,3). 

v  -  y 

The  two  polarization  directions  and  phase  velocities 
are  related  to  the  index  ellipsoid  in  the  following  manner. 
Assuming  an  arbitrary  direction  of  propagation  of  the  inci- 
dent  light  wave,  specified  by  the  unit  vector  s,  the  plane 

A  _ 

normal  to  s,  containing  D,  will  intersect  an  index  ellip¬ 
soid  of  the  medium  in  an  ellipse  as  shown  in  Fig.  2.  The 
radius  vector  r,  from  0  to  any  point  on  the  ellipse  can 
be  described  by, 

r 2  -  I  x2  (15) 

k«l  * 

where  each  x^  must  satisfy  Eg.  (11)  and,  using  Eq.  (10), 
the  relation 

D  •  s  *  |  x.  s.  -  0  ,  (16) 

k-1  *  K 

simultaneously.  At  four  points  along  the  ellipse  this  radius 
vector  will  be  an  extremum.  Obviously,  these  points  define 
the  axes  of  symmetry  of  the  section,  which  are  perpendicular 
to  each  other.  If  it  can  now  be  shown  that  when  Eq.  (5)  ap¬ 
plies,  r  is  an  extremum  of  the  index  ellipse  given  by 
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FIG  2  ELLIPSE  FORMED  BY  INTERSECTION  OF  PLANE  NORMAL  TO  S 

WITH  INDEX  ELLIPSOID  . 
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Eq.  (15),  then  the  relationship  between  the  light  polariza¬ 
tion  direction*  in  the  medium  can  be  clearly  established 
and  simply  determined  for  any  specific  case.  Furthermore, 
it  should  follow  from  the  nature  of  the  geometry,  that  the 
magnitude  of  r,  when  r  is  an  extremum,  will  be  propor¬ 
tional  to  the  appropriate  phase  velocity  vpi  or  vp£. 

The  necessary  condition  for  r  to  be  an  extremum  is 
obtained  by  *he  Lagrange  method  of  undetermined  multipliers. 

In  this  method  the  function 

* "  Ji  p  +  a“k*k  +  B(^T ' 11 

is  formed,  noting  that  the  last  two  terms  are  identically 
zero  by  Eqs.  (11)  and(l6).  Hence,  R  is  actually  eqrn  ta¬ 
lent  to  r2,  but  contains  the  constraints  on  r  in  the 
terms  having  the,  as  yet,  undetermined  multipliers  a  and 

0. 

If  the  derivative  of  R  is  now  taken  with  respect  to 
x^  and  set  equal  to  zero,  R  and,  therefore,  r  will  be 
an  extremum.  Performing  the  differentiation, 

*  (l  +  ^  GSfc  *  0#  (k  «  1,2,3)  (17) 

and  a  and  0  can  be  determined  by  the  following  simple 
considerations. 

Multiplying  Eq.  (17)  by  x*  and  summing  over  k  re¬ 
sults  in 


dR 
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3 

2 

k=l 


(1  +  +  aV’k 


3 

2 

k»l 


2  V 

X.  2  +  0~ 

*  S 


0, 


3 

since  2  -  0  by  Eq.  (16).  Applying  Eqs.  (11)  and 

k*l 

(15)  to  the  remaining  expression  then  yields. 


0  -  -r2. 


(18) 


On  the  other  hand,  multiplying  Eq.  v17)  by  s^  and 
again  summing  over  k  gives  immediately. 


3 

2 

k-1 

(1  +  +  “*k* 

3 

*  2 
k*l 

f-r2  Yk  + 

using  Eq. 

(18).  Bence, 

a 


3 

2 

k*l 


Vk 


s 


(19) 


Upon  substitution  of  Eqs.  (18)  and  (19)  into  Eq.  (17) 
the  condition  that  the  x^  must  satisfy  can  be  written  in 
terms  of  known  parameters  and  reduced  to  a  form  that  di¬ 
rectly  relates  to  the  optical  problem.  Thus, 


+  r2s.  2 
K  t=l 


xl*l 


0  (k 


1,^,3) 


is  the  condition  that  insures  that  r  is  an  extremum  of 
the  index  ellipse. 

To  transform  this  equation  into  electromagnetic  terms 
it  will  be  recalled  that 
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\  fA 

**'■<&  y[7  ’ 

by  Eqe.  (8)  and  (10).  Substituting  for  ^  then  gives 


.2  Dv  _  3 


^  ^  ■ o>  <k'l,2’3) 


(20) 


:k 


but 


D  •  D 


^  “  E  •  D 


by  Eqs.  (6),  (9),  (10)  and  (15).  Hencs,Eq.  (20)  becomes 


_2  D2s.  3 

a  -  _p  -)\  *  s,  v*  *  °- 

C  E  •  D  *  E  •  D  -t*l 


(21) 


Now,  final  reduction  may  be  effected  by  noting  that 
Eq.  (5)  can  be  written  as 


B-  -  E  -  (E  •  '»)»  -  (E  •  §)(§). 

n2 


with  reference  to  Fig.  1.  Therefore, 


n* 


<*  •  §)<§> 


E  •  D 


r2. 


(22) 


and  Eq.  (21)  simplifies  to 
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n  2  „  3 

(1  -  )°k  +  n  Sic  Z  "  °» 

k  l>*i 

or 

D*  -  na[Ek  -  (E  •  a)sk],  (k  -  1,2,3) 

which  is  just  the  component  form  of  Eq.  (5).  Thus,  when 
Eq.  (5)  applies,  the  light  waves  will  be  polarized  in  the 
directions  of  the  principal  axes  of  the  index  ellipse  and, 
therefore,  always  at  right  angles  to  each  other. 

The  relationship  between  the  magnitude  of  r  and  the 
phase  velocity  vp  is  established  by  Eqs.  (4)  and  (22), 
giving 


(23) 


Note  that  in  an  electrically  isotropic  medium,  where  the 
index  ellipse  degenerates  into  a  circle,  the  two  phase 
velocities  are  equal,  precluding  the  necessary  condition 
for  double  refraction. 

Summarizing  these  results,  it  is  seen  that  the  two 
monochromatic  plane  light  waves  that  propagate  in  an  elec¬ 
trically  anisotropic  medium  are  linearly  polarized  at  right 
angles  to  each  other,  with  their  directions  of  polarization 
coincident  with  the  principal  axes  of  the  index  ellipse 
formed  by  the  intersection  of  the  plane  normal  to  the  direc¬ 
tion  of  light  propagation  and  the  index  ellipsoid  of  the 
medium.  The  phase  velocities  of  the  waves  are  inversely 
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proportional  to  the  half-lengths  of  the  principal  axes  of 
the  index  ellipse.  When  che  half-lengths  are  equal ,  the 
net  birefringence  is  zero.  This  condition  occurs  only  when 
the  light  is  propagating  either  along  the  optic  axis  of  a 
crystal  or  in  an  electrically  isotropic  medium. 

B.  elastic  properties  of  amorphous  solids 

As  previously  indicated,  an  amorphous  solid  must  be 
stressed  in  order  to  create  the  conditions  for  birefringence. 
The  primary  mechanisms,  however,  that  transform  the  internal 
structure  of  the  material  from  the  isotropic  to  the  anisotro¬ 
pic  state  are  the  elastic  deformations  or  strains  induced  by 
the  stress. 16  It  is  therefore  natural,  then,  to  develop  the 
subsequent  elasto-optic  interaction  theory  in  these  terms, 
rather  than  in  terms  of  the  stresses,  in  order  to  provide  a 
true  picture  of  the  phenomena.  Once  the  formulation  is  com¬ 
plete,  the  conversion  to  the  more  useful  stress  form  is 
easily  accomplished  by  means  of  the  stress-strain  relation¬ 
ship. 


6 

s 

-i=i 


CklSl 


(k  ■  1  to  6) 


(24) 


where  T^  are  the  components  of  stress,  the  components 

of  strain  and  the  elastic  constants  of  the  medium. 

The  complete  symmetry  of  the  amorphous  solid's  internal 
structure  reduces  the  number  of  independent  elastic  constants 
to  two,  C  and  C  However,  a  shear  '>ntant,  C  , 

A A  44 

exists,  which  may  be  expressed  in  terms  of  t.ie  independent 
constants  as. 


C 

44 


c  -  c 
—11 _ 

2 


(25) 
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The  matrix  of  the  elastic  constants  for  this  type  of  material 
is  then 


with  no  coupling  between  normal  and  shear  components. 


Consider  now  an  initially  strain-free  amorphous  solid 
in  the  form  of  a  plate  whose  major  faces  are  parallel  to 
the  x^x^-plane,  as  shown  in  Fig.  3»  The  light  wave  pro¬ 
pagates  in  the  -x  -direction  and  the  elastic  wave  is  con- 

3 

strained  to  propagate  in  the  x^x^plane  only.  The  dimen¬ 
sions  of  the  plate  are  assumed  to  be  large  enough,  compared 
to  the  elastic  wavelength,  so  that  the  wave  propagates 
freely  in  the  medium,  i.e.,  uninfluenced  by  the  boundaries. 

A 

Let  k  be  the  unit  vector  specifying  the  direction  of 
elastic  wave  propagation  in  the  x^x^-plane,  and  let  <t> 
be  the  angle  this  vector  makes  with  the  x^-axis. 

If  the  elastic  wave  is  a  plane,  harmonic  s'ave  of  fre¬ 
quency  fe,  then  the  elastic  displacement  at  a  point  P 
in  the  plate  is  ^ iven  by 


u  =  U  cos(o>et  -  k  •  r) 


(27) 


T-3/321 


-17- 


COLUMBIA  UNIVERSITY— ELECTRONICS  RESEARCH  LABORATORIES 


where  U  is  the  displacement  amplitude,  cd^  *•  2 k  «  k 
is  the  propagation  vector,  is  the  elastic  wavelength, e 

and  r  is  the  position  vector  to  the  point  P. 

In  order  for  Eq.  (27)  to  represent  a  true  displacement 
in  the  elastic  medium, it  must  satisfy  the  equations  of 
motion,  which  in  their  most  general  form  are, 17 


P 


yuk 

dt= 


3  d2ut 
t-1  ktd*2 


(k 


1,2,3)  . 


(28) 


Here,  u]c  are  the  components  of  u,  p  the  density  of  the 

medium,  6^  elastic  moduli  that  are  functions  of 

and  <*>,  and  £  a  new  variable  related  to  k  •  r  through 


k  •  r 


(29) 


The  elastic  moduli,  6^,  for  an  amorphous  solid 
supporting  a  wave  propagating  in  the  x^x^-plane  only,  are* 


6 

ss 

a2C 

+  02C 

li 

n 

44 

6 

6  ■ 

a8(C  +  C  ) 

12 

21 

12  44 

6 

* 

6  = 

0 

13 

31 

6 

- 

a2C 

+  02C 

22 

44 

11 

6 

ar 

6  = 

0 

23 

32 

6 

SB 

C 

33 

44 
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where  a  *•  cos  <$>,  P  -  sin  <t>  are  the  direction  cosines  of 

with  respect  to  the  x  ,x  axes. 

1  r 

Substituting  these  moduli  into  Eq.  (28)  defines  the 
specific  equations  of  motion  relevant  to  the  geometry  of 

Fig.  3# 


d2U 


dt2 

d2u, 

dt2 

d2u 


Jl  . 


d2u 

d2U 

(a2C 

+  d2C  ) - 1 

+  a0(c 

+  c  ) — -f- 

n 

44  de2 

12 

44  d*2 

d2U 

d2u 

a0(c 

+  c  ) - 1  + 

(a2C  + 

d2c  ) — * 

'  12 

44  d*2 

44 

11  d^2 

(30) 


_  d2^ 

P  dt2  *  44  d$2 


If  Eq.  (27)  is  now  written  in  component  form  as. 


uJt  -  Uj^  cos(o>et  -  ~£)  ,  (Jc  -  1,2,3) 


using  Eq.  (29),  and  the  operations  indicated  in  Eq.  (30) 
are  performed,  there  results  the  following  se*  of  conditions 
on  the  U,  : 

r(a2c  +  e2c  )  -  pv_2iu  +  ad(c  +  c  )u  »  o 

1 V  11  *4  e  J  1  12  44  2 

ae(Ci2  +  C44)01  +  t(o*C44  ♦  82Cit)  -  pve2]U2  -  0  (31) 

(C«  ‘  pveS)U3  ■  0 

where  v  is  the  speed  of  propagation  of  the  elastic  wave, 
e 
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These  conditions  allow  the  determination  of  the  speed 
of  propagation  and  the  relationships  between  the  U^.  In 
the  simplest  case,  letting  U  ■  U  ■  0  and  U  ^  0,  the 
last  of  Eg.  (31)  yields. 


Since  the  direction  of  propagation  is  in  the  x^x^-plane, 

U  is  transverse  to  this  direction  and  thus  represents  the 

3 

propagation  of  a  pure  transverse  wave  whose  displacements 
are  parallel  to  the  direction  of  light  propagation,  s. 

When  U  ,  U  /  0  and  U  ■  0,  the  displacements  are 
i  2  3 

confined  to  the  x^x^-plane  and  the  relationship  between 
them  can  be  obtained  by  noting  that  unique  values  of  v 

Q 

will  exist  only  when  the  determinant  of  the  coefficients 
is  equal  to  zero.  Setting  the  determinant  equal  to  zero, 
and  solving  for  v  ,  gives  two  distinct  solutions  to  the 
determinants!  equation. 


By  substituting  these  expressions  back  into  the  first  of 
Eq.  (31)  it  is  found  that, 

U  =  *7  U 
2  CL  1 
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with  v  •  v  .  Therefore,  by  a  simple  geometric  interpre- 
e  e 

2 

tat ion  of  these  relationships,  it  is  established  that 

x 

is  the  spesd  of  propagation  of  a  longitudinal  wave,  while 

v  is  tM  propagation  speed  for  a  transverse  wave  whose 
2 

displacements  are  always  perpendicular  to  those  of  the  longi¬ 
tudinal  wave.  Hence, 

(32) 

and  it  may  be  concluded  that  an  amorphous  solid  will  support 
two  types  of  waves  simultaneously,  a  pure  L-wave  with  dis¬ 
placements  in  the  direction  of  k  and  a  pure  T-wave  with 
displacements  normal  to  k. 

C.  THE  ELASTO-OPTIC  INTERACTION 

In  the  strain-free  plate  of  Fig.  3>  as  already  noted, 
the  index  ellipsoid,  Eq.  (11),  at  any  point  P  is  degen¬ 
erate  to  a  sphere  because  the  medium  is  optically  (i.e., 
electrically)  isotropic.  The  radius  of  this  sphere  is 
equal  to  where  e  is  the  characteristic  dielectric 

constant.  Hence,  by  Eqs.  (11)  and  (22),  the  expression 
for  the  circular  section  in  the  x  x  -plane  may  be  written 

X  2 

in  the  form. 


e 


1. 


(j3) 
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Wien  an  elastic  wave  is  propagating  in  the  plate,  the 
induced  strains  modify  the  dielectric  properties  of  the 
medium  so  that  it  becomes  birefringent.  This  birefringent 
state  is  expressed  geometrically  by  distorting  the  circular 
section  into  an  ellipse  which,  in  general,  has  its  principal 
axes  rotated  through  an  angle  0,  as  shown  in  Fig.  4.  The 
equation  of  the  index  ellipse  is 


a  x2  +  a  x2  +  2a  xx  *  1  , 

U  1  222  12  1  2 


(34) 


where  the  are  known  as  polarization  constants.  They 

will  be  shown  later  to  be  reciprocal  to  the  squares  of  the 
refractive  indices  in  the  index  ellipse. 

The  fundamental  assumption  of  the  linear  theory  of 
birefringence  in  solids  is  that  the  difference  between  the 
polarization  constants  of  Eq.  (34)  and  the  corresponding 
reciprocals  of  the  square  of  the  characteristic  refractive 
index  of  Eq.  (33)  i*  a  linear  combination  of  the  components 
of  strain.  In  general,  there  will  be  six  equations  in  six 
variables.  However,  because  the  constants  in  the  linear 
combination  display  the  same  matrix  pattern  as  the  elastic 
constants  in  Eq.  (£6),  the  number  of  terms  in  these  elasto- 
optic  equations  will  reduce  to: 


a  -  — -  ■  p  S  +  p  S  +  p  S 

11  n2  11  11  12  22  ’2  33 

0 


a  _  apS  +pS  +pS 
22  n2  12  11  11  22  12  33 

0 


a 

12 


P  S 

44 


12 


(35) 
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FIG.  4  PROPERTIES  OF  INDEX  ELLIPSES  IN  STRAINED  AND  UNSTRAINED 
STATES  FOR  ELASTIC  WAVE  PROPAGATION  IN  X,-  DIRECTION  AND 
3/32!  UGHT  WAVE  PROPAGATION  IN  X3- DIRECTION  .  _2if_ 
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where  consideration  has  been  given  to  the  fact  that  varia¬ 
tions  in  the  index  ellipsoid  in  planes  normal  to  the  x1x2“ 
plane,  for  the  geometry  of  Fig.  4,  have  no  effect  upon  the 
birefringence. 

The  in  Eq.  (35)  are  the  Pockel  e  las  to-  optic  con¬ 

stants  for  an  amorphous  solid  and  the  Sj^  an e  the  strain 
components.  The  elasto-optic  constant  p  is  not  an  in¬ 

dependent  constant,  but  is  related  to  the  other  two  by  an 
expression  similar  to  Eq.  (25), 


(36) 


This  is  due  to  the  effects  of  symmetry,  where  it  can  be 
shown  that  in  an  amorphous  medium  only  two  independent  con¬ 
stants  are  necessary  to  describe  the  elasto-optic  behavior. 

The  components  of  strain  in  Eq.  (35)  are  functions  of 
the  displacement  components  according  to  the  following 
formulas  from  elementary  elasticity  theory: 


du 


(i  /  <), 


(37) 


with  the  x^  here  as  spatial  coordinates.  In  these  terms, 
using  Eq.  (27)* 


» 


I 
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S  »  T2  U  cos  <*>  sin  (cot  -  k  •  r) 
a  *e  1  e 


S  -  ~  u  sin  0  sin  (co  t  -  k  •  r) 

22  A.  2  e 

(38) 


S  =  0 
33 


S  m  ~  (u  sin  4>  +  U  cos  <f)sin(aj  t  -  k  •  r)  . 
12  Afi  1  2  e 


Equation  (38)  may  be  put  into  a  more  useful  form  by 
converting  it  to  functions  of  and  U^,,  i.e.,  longi¬ 

tudinal  and  transverse  wave  displacements  (Fig.  5) •  There 
results  t 


L- Waves 


S  ■  T~~  U_  cos z<t>  sin(u>  t  -  k  •  r) 

11  XeL  L  e 

S  «  — -  UT  sin z<t>  sin(co_  t  -  k  •  r)  (39) 

22  XeL  t  e 

S  *  --  U.  sin  2<f>  sin(co_t  -  k  *  r)  . 

12  XeL  L  e 

T- Wave  s 

S  *  U_  sin  <t>  cos  0  sin(co  t  -  k  •  r) 

11  \.T  T  e 

s  ■  u_  sin  <t>  cos  <t>  sin(co  t  -  k  •  r)  (40) 

22  AeT  T  e 

S  ■  Um  cos  2 <t>  sin(o)  t  -  k  •  r)  . 

12  A„m  T  e 
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FIG. 5  GEOMETRY  FOR  TRANSFORMATION  OF  ELASTIC 
DISPLACEMENT  COMPONENTS 
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With  the  establishment  of  the  elasto-optic  equations 
and  the  strains  that  are  pertinent  to  Fig.  4,  it  is  now 
possible  to  proceed  directly  to  a  determination  of  the  angle 
6  and  the  refractive  index  changes  along  the  principal  axos 
of  the  index  ellipse. 

The  angle  0  is  obtained  by  making  a  coordinate  trans¬ 
formation  from  the  x^x^  axeB  to  the  xi,xi  axes  which 
are  coincident  with  the  principal  axes.  The  transformation 

equations 


x»  ■  x  cos  ©  +  x  sin  0 
11  2 


x'  «-x  sin  ©  +  x  cos  9 
2  1  2 


are  first  substituted  into  the  equation 


ax,2+ax,2*a  x2  +  a  x2  +  2a  xx 
11  22  11  1  22  2  12  1  2 


which  describes  the  invariance  of  the  index  ellipse  to  such 
a  transformation.  Then,  the  coefficients  of  like  variables 
are  equated  to  yield. 


a  *  a  cos2©  +  a  sin2© 
111  2 


22 


m  a  sin2©  +  a  cos2© 


(41) 


2a  -  (a  -  a  )  sin  2©  . 
12  1  2 


If  the  second  of  Eq.  (41)  is  now  subtracted  from  the  first, 
and  the  third  divided  by  sin  20, 


T- 3/321 


-28- 


COLUMBIA  UNIVERSITY— ELECTRONICS  RESEARCH  LABORATORIES 


a 

1 


a 

2 


a  -  a  2a 

— Li _ Z2.  ^  12 

cos  20  sin  2 9  ’ 


Hence ,  the  orientation  of  the  index  ellipse  is  governed  by, 

2a 

tan  20  =  - - ^ —  .  (42) 

11  22 

Substitution  of  Eqs.  (39)  and  (40)  into  Eq.  (35)  leads 
to  specific  forms  of  the  polarization  constants  in  Eq.  (42). 
Then,  9  may  be  evaluated  for  both  L  and  T  waves.  It  is 
found  that. 


(43) 


and 


tan  20t  •  -  cot  2 <t>  . 

The  latter  expression  can  be  put  into  a  simpler  form 
by  writing  it  as, 

cos  20t  cos  2 <t>  +  sin  20T  sin  2 <t>  -  cos  2(0T  -<*>)«  0, 
and  recognizing  that  in  order  to  satisfy  this  equation, 

2(eT  -  0)  -  §  . 

Thus, 

0T  “  ♦  +  f  •  (44) 
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The  L  wave  is  therefore  seen  to  polarize  the  light 
in  the  k  direction,  and  at  right  angles  to  this  direction. 
A  T  wave,  on  the  other  hand,  splits  the  ircident  light 
into  two  components  that  make  angles  of  45  and  135  deg  with 
the  k  direction.  Since  Eqs.  (43)  and  (44)  are  functions 
of  <t>  only,  the  orientation  of  the  index  ellipse  is  also 
seen  to  be  the  same  at  every  noint  in  the  elastic  wave  field. 

Although  the  orientation  of  the  index  ellipse  is  fixed 
once  the  direction  of  elastic  wave  propagation  is  known,  the 
lengths  of  the  principal  axes  will  vary  with  the  variations 
in  strain.  Let  the  change  in  length  of  these  axes  be  de¬ 
ne,  ted  by. 


fir 


l 


(45) 


as  indicated  in  Fig.  4.  Since  the  phase  shifts  in  the 
transmitted  light  are  of  major  interest,  it  is  important  to 
express  fir  and  fir^  as  changes  in  refractive  index,  in 
terms  of  the  elasto-optic  properties  of  the  medium  and  the 
strains. 

Toward  this  end,  it  is  first  noted  that, 

a  +  a  *  a  +  a 
1  2  11  22 

2a 

ai  as  sin  29  * 
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from  Eq.  (4l).  Hence, 


a 

l 


a  -fa 

J.1  ,  22 

2 


a 

*  sin  29 


a  +  a 

_1 2 _ 22 

2 


a 

_ 12 _ 

sin  29  * 


(46) 


If  Eq.  (35)  is  now  substituted  for  the  polarization  con¬ 
stants  and  the  resulting  expressions  applied  to  Eq.  (45) , 


l 


(p  +  p  )  (S  +  S  ) 
-  JU. _ la  11  12 


+ 


P 


S 

4  .12 


8  in  29 


(P 


UL_ 


p  )  (n 

12  .  JU 


Sig) 


p  s 
lai-li 

sin  2f 


These  formulas  can  be  simplified  by  forming  the  sum 
and  difference  functions. 


6r  +  6r 


(P 


li 


P  )  (s 
12  11 


S  ) 
12 


6r 


l 


2p  S 

_ 44  12 

sin  29  > 


(47) 


and  eliminating  the  sin  29  as  follows. 

When  Eq.  (42)  is  written  in  terms  of  the  strain",  by 
substitution  of  Eq.  (35) > 


tan  29 


2p  S 

_ ,14  12 . 

-  p  )(S 
lit  '  11 
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This  equation  has  a  geometrical  interpretation.  The 
numerator  is  just  the  altitude  of  a  right  triangle  whose 
hypotenuse  makes  an  angle  29  with  the  base,  or  denomina¬ 
tor.  The  length  of  the  hypotenuse  is,  therefore. 


h 


r(S 


11 


-  s  )2  +  4p  2s  2  . 

22  44  12 


Then,  the  sin  29  mu<*t  be. 


2p  S 

sin  29  -  — , 


and  the  difference 


fir  -  fir 


(P  -  P  ) [ (S  -  S  )2  +  S  2] 
11  12  11  22  12  1 


(48) 


using  Eq.  (36).* 

A 

Considering  first  an  L  wave  propagating  in  the  k 
direction,  Eqs.  (47)  and  (48;  now  allow  a  unique  representa¬ 
tion  of  fir^  and  fir^,  in  their  simplest  form,  once  the 
stains  are  written  in  terms  of  the  displacement  UL  by 
means  of  Eq.  (39)»  The  sum  and  difference  functions  in  this 
case  are. 


fir  +  fir 
1  2 


(p  +  P  )U. 
11  12  L 


(49) 


fir  -  fir 
1  2 


(p  -  p  )UT  , 
11  12  I* 
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neglecting  the  common  multiplicative  factor.  If  P1X  >  P12 

and  S  is  of  positive  sign,*  the  solutions  to  Eq.  (49) 

12 

are : 


fir 

1 


P  UT 
11  L 


6r 

2 


P  UT 
12  U 


(50) 


With  S  negative,  an  apparent  rotation  of  che  principal 
12 

axes  through  90  deg  occurs  and. 


fir 


P  Ur 
12  L 


fir 


P  Ut 
11  I* 


The  solutions  for  a  T  wave  follow  just  as  easily  by 
application  of  Eq.  (40).  Here,  the  sum  and  difference  func¬ 
tions  are. 


fir  +  fir  =  0 
1  2 


fir  -  fir 


+  (p 
-  11 


PjUT 


Therefore, 


fir  =>  -fir  *  +  p  U_  .  (51) 

1  2  ~  44  * 

The  most  useful  forms  of  Eqs.  (50)  and  (51)  ar®  obtained 
when  <t>  *  0,  i.e.,  the  elastic  wave  propagates  in  the 

2p..*s,*> 

*  The  sign  of  S  is  important  because  fir^  -  fir^  -  si*4202 
s  h  is  positive  or  negative  depending  on  the  sign  of 
when  0  <  6  <  n/2  and  p++  >  0. 
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x  -direction.  Then,  Eq.  (39)  indicates  that  only  /  0, 

and  Eq.  (40),  that  only  t  0.  Recalling  the  multiplica¬ 

tive  factors  implicit  in  Eqs.  (50)  and  (51) >  it  is  seen  that 
these  formulas  now  become: 

L- Waves 


6r  *  p  S 

i  *ii  li 


6r  =  p  S 

2  12  11 


(52) 


T- Waves 


6r 


i 


P  S 

44  12 


(53) 


Eqs.  (52)  and  (53)  are  working  equations  which  will  be  used 
exclusively  in  what  follows. 

The  refractive  index  changes,  corresponding  to  the 
changes  in  length  of  the  principal  axes  indicated  in  Eqs. 

(52)  and  (53),  may  be  derived  by  noting  the  relationship 
between  n  and  r  in  Eq.  (22)  and  the  forms  of  Eqs.  (33) 
and  (34).  It  will  be  recalled  that  in  the  development  of 
Eq.  (22),  the  index  ellipsoid  was  given  by  Eq.  (11)  which 
is  of  the  same  form  as  Eq.  (33).  Therefore,  it  is  concluded 

that 


1  1 
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1  2  p 

n  n  *- 
1  o 


(54) 


2  n  2  n  2 


by  Eq.  (45). 


Applying  Eq.  (54)  to  Eq.  (52)  now  gives. 


n  -  n  \  ,n  +  n 


n  2n  2 


»n  -  n 


■■  1 1  *  p  S 

3  11  11 


O  1 


-n  -  n  \  ,n  +  n 


n  2n  2 


in  -  n 


P  s 
12  11 


O  2 


if  n^  and  n^  differ  only  slightly  from  n  .  The  elasto- 
optic  constants  on  the  right  are  usually  expressed  in  terms 
of  the  so-called  Neumann  strain-optic  constants  for  isotro¬ 
pic  media. 


li  n 


12  n 


(55) 


hence  the  above  formulas  may  finally  be  written  as. 


6n  .  =  -n  2qS 

o  ^  11 


fin  T  =  -n  2pS 
2"  0  11 


(56) 


where  fin  .  =  n  -  n  . 

liSL*  1  >  2  0 
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The  difference  in  the  refractive  index  changes  in  Eq.  (52), 


1  1  «.  2(ns  "  ni)  /_  _  No 

*ri  -  6r2  "  - - ;  ~  \  3  (pn  Pxe,Sn 

1  z  n  z  n  nJ 

12  o 


defines  the  strain-optic  constant  b  given  in  the  literature. 


6n  n  -  n 

» s  — 1  -  n0  (q  * p) 
ii  ii 


(57) 


where  fin  2  n  -  n  and  Eq.  (55)  is  used.  This  constant 
21  2  l 

simplifies  the  expressions  for  T  waves,  corresponding  to 
those  in  Eq.  (56). 


The  development  for  T  waves  proceeds  in  similar  fash 
ion  to  that  used  for  L  waves  except  that  Eq.  (53)  re¬ 
places  Eq.  (52).  Thus, 


n  -  n 
i  o 


P44S 


12* 


Now,  introducing  Eqs.  (36)  and  (55)  into  this  equation  puts 
it  into  the  form. 


n  -  n  ■  -(n 
10  2 


nft2(q  -  p) 

_q  _ 


S 

12 


12* 


Hence , 


fin 


iT 


(58) 
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The  conversion  of  Eqs.  (56)  end  (58)  to  the  stress  forms, 
which  will  find  immediate  application  in  the  solid  light  modu¬ 
lators,  is  accomplished  by  first  inverting  Eq.  (24)  and  ex¬ 
pressing  the  coefficients  in  terms  of  the  commonly  known  pro¬ 
perties  of  an  isotropic  medium.  It  can  be  shown  that,18 

>  <*•<•  -  1.2.3) 

exactly  when  k  i-  l  and  approximately  when  k  ■  <t.  For 
practical  purposes  this  equivalency  is  satisfactory  in  nil 
cases,  with  o  the  Poisson's  ratio  and  E  the  modulus  of 
elasticity  of  the  medium. 

Upon  substitution  for  the  in  the  elasto-optie 

equations  the  changes  in  refractive  index  are  then  deter¬ 
mined  by  the  formulas: 


6nlL  - 


O  il 


6\l  *  -<1_r!I>no£pTii  • 


(59) 


and 


SV  *  *Sn2T  * 


-It 

2  12 


(60) 


where  B,  the  stress-optic  constant,  is  defined  by 


B  =  (1  p  ?)b  »  (-1  |  -g)noa(q  -  p) . 


(61) 
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III. 


ION  BY 


AMORPHOUS  SOLIDS, 


As  a  result  of  a  stress  wave  propagating  in  an  amorphous 
solid  it  has  been  shown  that  refractive  index  changes  wi.l 
occur  in  specific  planes  of  polarization,  at  each  point  in 
the  stress  wave  field,  to  an  extent  dependent  upon  the  mag¬ 
nitude  and  direction  of  the  stress  at  that  point. 

Obviously,  the  refractive  index  changes  alone-  the  direc¬ 
tion  of  elastic  wave  propagation  will  conform  to  the  shape 
of  the  traveling  wave,  as  is  evident  from  Eqs.  (59)  and  (60). 
Thus,  in  effect,  a  variable  optical  transmission  gracing  is 
created  whose  configuration  is  directly  related  to  the  char¬ 
acter  of  the  stress  wave  at  each  instant  of  time. 

If  the  variations  in  refractive  ind.*x  are  small  and  the 
thickness  of  the  medium,  in  the  direction  of  light  propaga¬ 
tion,  is  such  that  internal  bending  (internal  refraction) 
can  be  neglected,*  the  diffraction  effects  caused  by  sinus¬ 
oidal  excitation  of  the  solid  c  an  be  explained  by  the  simple 
Raman-Nath  theory14. 

consider  the  geometry  in  Fig.  6.  The  light  wave  propa 
gates  along  the  x3-axis  in  the  direction  specified  by  s. 
Its  i  vector  is  polarized  at  an  angle  7  with  respect  to 

the  XjL-axif , 

A  plane,  harmonic  elastic  wave  is  assumed  to  propagate 
Jn  the  XjX, -plane,  along  the  x^axls,  in  the  direction 
of  f..  It  induces  an  index  ellipse  whose  major  axis  makes, 
in  general,  an  angle  0  with  the  Xj-axis. 


•  The  effect  of  internal  refraction  on  the  diffraction  pat- 
tern  is  discussed  in  tiec,  IV-D, 
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l-IG.  6  DIFFRACTION  GEOMETRY 
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When  the  light  wave  enters  the  stressed  solid,  it  splits 
into  two  waves  which  undergo  phase  shifts  with  respect  to  the 
phase  of  the  incident  wave  according  to. 


*1 


2tt 


where  is  the  wavelength  of  the  light  and  6d  is  the 

change  in  optical  path  length.  The  optical  path  length 
changes  are  related  to  the  changes  in  the  refractive  indices 
by  the  equation, 

6d  =  L6n 

vhere  L  is  the  geometrical  path  traversed  in  the  elastic 
wave  field. 

Since  the  elastic  wave  is  harmonic,  the  amplitudes  of 
the  phase  shifts  are  of  primary  concern.  This  amplitude  is 
known  as  the  modulation  index  and  is  defined  by. 


2r 

^  5n  L, 

rm  \  m  9 


(62) 


with  6n  the  maximum  change  in  the  refractive  index  and 
m 

L  now  interpreted  as  the  width  of  the  elastic  wave  fi^ld 
in  the  direction  of  s. 

Substituting  Eq,  (59)  into  Eq,  (62),  the  modulation  in¬ 
dices  for  an  L  wave  are  obtained t 


K  -  rI(HJi>no“'IL  T, 


t 


m 


iL  \ 

-  n„*pL  T, 


mL 


(63) 


2L 


V 


mL 
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where  T  ,  is  the  maximum  normal  Strega  T  • 
ml  A 


For  a  T  wave, 

t 


m 


liT 


V 


m 


lj>T 


(64) 


upon  substitution  of  Eq.  (60),  with  maximum  T1#  denoted  by 

T  • 

ms 

Consider  now  that  the  plate  of  Fig.  6  is  uniformly  il¬ 
luminated  by  a  perfectly  collimated,  linearly  polarized  beam 
of  light.  This  beam  will  be  modulated  by  the  elasto-optic 
interaction  in  the  plate  so  that,  at  the  output  face,  its 
wavefront  will  vary  sinusoidally  along  the  x^axis.  If  a 
lens  of  focal  length  f  is  interposed  between  the  plate  and 
a  screen  which  is  placed  in  the  focal  plane  of  the  lens,  a 
diffraction  pattern  in  the  form  of  spots  of  light  will  be 
observed  on  this  screen.  The  position,  with  respect  to  the 
optical  axis,  of  these  spots  or  diffraction  orders,  along  a 
line  parallel  to  £,  is  given  by  the  formula,* 


<Lj i  B  f  tan  fi^,  (h  *=  0,1,2,,,.) 

where  the  diffraction  angle  0^  is  defined  by 


V 

sin  ^  ~  . 


When  \  «  *  ,  sin  ^  ^  in  the  lower  orders,  and 


**  -  f  “ k  f 


e 


(65) 


is  sufficiei  tly  accurate  for  moat  computational  purposes. 


*  The  diffraction  pattern  in  the  direction  normal  to  £  will 
be  neglected  because  it  is  not  pertinent  to  these  considera¬ 
tions  , 
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The  light  intensity  distribution  in  the  focal  plane  of 
the  lens  will,  in  general,  depend  on  the  angle  y  of  J:he 
E  vector,  ae  well  as  the  modulation  indices*  When  E  is 
polarized  in  the  direction  of  one  of  the  principal  axes  of 
the  index  ellipse,  Raman-Nath  theory  indicates  that  the  max¬ 
imum  value  of  the  electric  intensity  in  the  k**1  order  will  be 

^  (66) 

where  is  the  Xth  order  Bessel  function.  Hence,  the  max¬ 

imum  relative  light  intensity  in  the  Xth  diffraction  order  is* 

No  change  in  the  light  polarization  will  occur  in  this  case 
because  of  the  initial  alignment  with  the  principal  axes. 

In  the  more  general  case  where  the  E  vector  xs  re¬ 
solved  into  birefringence  components, 

E  ■  E  cos(0  -  y) 

1  (67) 

z  -  E  sin(S  -  y), 
s 

the  situation  becomes  more  complex. 

First  of  all,  each  of  the  above  components  gives  rise  to 
a  diffraction  pattern  where  the  maximum  electric  intensities 

are, 

*  The  relative  light  intensity  distribution  in  the  k^1  order 
is  of  the  form  sin2j0^«  since  sinc  x  mul¬ 

tiplier  is  of  no  importance  to  the  present  discussion,  all 
formulas  for  relative  intensity  will  be  concerned  with  the 
maximum  relative  intensity  only.  A  complete  discussion  of 
the  relative  light  intensity  distribution  to  be  expected 
in  the  above  situation  can  be  found  in  Ref.  12. 
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S  v  ■  E  J.  (0  ) 

ik  i  kvrmi/ 

E  -  E  J.  (0  ). 

2k  2  kvrm ZJ 


(68) 


These  intensities  will  be  deviated  through  the  same  diffrac¬ 
tion  angle  and  since  they  originate  from  the  sane  in¬ 

cident  light,  by  diffraction  on  the  same  elastic  wave,  the 
resultant  electric  intensity  in  the  kth  order  will  then  be 
Just  their  vector  sum.  Thus,  linear  polarization  is  pre¬ 
served;  however,  the  polarization  of  the  diffracted  light  can 
now  differ  from  that  of  the  incident  light. 


The  angle  of  polarization  of  the  diffracted  light  is 
easily  determined  by  considering  the  geometry  in  Fig,  7  and 
Eg.  (68).  If  the  ratio 


E„k  E  J.  (0  ) 

is  formed  and  Eq.  (67)  applied,  the  polarization  angle  is 
given  by 

t,n(9  -  7k)  “  tan(9  -  7)  . 

For  L  waves,  recalling  Eq.  (43)  with  0-0, 

tan  7k  -  tan  7  tS/S  i  (69) 

k'nm ' 

while  for  T  waves,  the  angle  becomes, 

tan(f  "  7k)  “  (-Dk  tan(jf  -  y),  (70) 
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once  it  is  recognized  that 

Jk('lWr)  "  Jk<-*nuT>  ■  (-1)k  T)  . 
and  Eq,  (44)  is  employed. 

The  maximum  relative  light  intensity  can  also  be  derived 
immediately  by  referring  to  Fig,  7  where 

g  2  £  2 

=  E“  *k  =  Vty mi >  coe‘(e  -  y)  +  8i"2 -  7> 

upon  substitution  of  Eqs.  (67)  ard  (68).  In  the  two  cases  of 
interest,  this  equation  becomes: 

W  -  V<W  coe‘y  +  s1"27  (71) 

XkmT  “  Jk‘<*m1T)  (72) 

where  it  is  seen  that  for  T  waves,  the  maximum  relative 
light  intensity  in  all  orders  is  independent  of  the  polar¬ 
ization  of  the  incident;  wave. 
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IV,  APPLICATION  TO  SOLID  LIGHT  MODULATORS 


The  results  of  the  theory  of  birefringence  in  amorphous 
solids  will  now  be  applied  to  design  considerations  for  de¬ 
velopment  of  a  wide-band  (50  MHz)  solid  light  modulator  op¬ 
erating  at  a  center  frequency  of  100  MHz.  These  considera¬ 
tions  represent  the  initial  effort  in  the  solid  light  modu¬ 
lator  research  program  that  is  presently  being  conducted  at 
CUERL. 

The  modulator  is  essentially  an  ultrasonic  delay  line 
in  the  form  of  a  plate  of  very  high  optical  and  acoustical 
quality  fused-silica.  At  one  end  of  the  plate,  elastic 
waves  are  generated  by  quartz  transducers  resonating  in  the 
fundamental  mode.  At  the  other  end,  a  highly  efficient  ab¬ 
sorber  attenuates  the  waves  to  a  level  (  >  40  dB)  where  the 
delay  line  appears,  for  all  practical  purposes,  to  have  an 
infinite  insertion  loss.  Thus,  the  elastic  waves  may  be 
considered,  in  an  analytical  sense,  as  pure  traveling  waves. 

Illumination  of  the  modulator  in  the  experimentation 
will  be  provided  by  a  high- intensity  gas  laser  source  which 
emits  a  beam  of  light  with  a  0.6328^  wavelength.  The  light 
is  highly  collimated  over  the  optical  aperture  by  a  preci¬ 
sion  optical  system  so  that  a  plane  wave  normally  incident 
on  the  plate  is  achieved. 10 

Although  the  objective  of  the  program,  as  stated,  is 
the  development  of  a  100  MHz  device,  the  state-of-the-art 
in  fabrication  of  very  high  frequency  fundamental  mode  quartz 
transducers,  which  meet  the  program  requirements,  is  such, 
that  this  frequency  appears  to  be  at  the  outer  limit  of  the 
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fabricator's  capability.  Hence,  in  what  follows,  considera¬ 
tion  will  also  be  given  to  80  MHz  modulators  since  a  higher 
confidence  level  exists  in  the  fabrication  of  delay  lines 
at  this  frequency. 

A*  SOLID  LIGHT  MODULATOR  CONFIGURATIONS 

The  simplest  solid  light  modulator  is  shown  in  Fig.  8. 
On  the  left  end  of  the  delay  line  is  the  quartz  transducer 
crystal,  attached  to  the  line  under  high  pressure  and  tem¬ 
perature  by  either  ar.  indium  or  lead  bond.  This  "bond”  also 
serves  as  the  ground  electrode  for  the  transducer  electrical 
connections.  The  positive  electrode  covers  the  exposed  face 
of  the  crystal  and  determines  the  transducer  size,  for  only 
that  area  of  the  crystal  between  the  two  electrodes  will  be 
exc^^e<^  piezoelectrically,  This  electrode  may  be  a  tit  in 
vacuum  deposited  film  of  indium  or  a  conductive  foil  sup¬ 
ported  on  a  viscous  fluid  adhesive  layer,  like  vaseline  or 
vacuum  grease. 

At  the  right  end,  a  wedge-shaped  acoustical  absorber 
is  cut  into  the  line.  The  angle  of  the  wedge  depends  upon 
the  required  amount  of  attenuation  and  the  achievable  energy 
loss  per  reflection  when  the  edges  of  the  wedge  are  coated 
with  lead,  a  material  closely  matching  the  mechanical  im¬ 
pedance  of  f used-silica  for  shear  waves. 

By  a  geometrical  analysis,  this  angle  is  given  by 


CD 


-180° 
n  +  1  » 


where  n  is  the  total  number  of  reflections  in  the  wedge. 
If  n  is  even,  the  analysis  indicates  that  the  reflected 
beam  will  emerge  from  the  absorber  parallel  to  the  inclined 


T- 3/321 


-47- 


COLUMBIA  UNIVERSITY-ELECTRONICS  RESEARCH  LABORATORIES 


FIG  8  SINGLE  TRANSDUCER  SOLID  LIGHT  MODULATOR  CONFIGURATION 
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edge.  If  n  is  odd,  the  reflected  beam  will  emerge  paral- 
lei  to  the  incident  beam.  Thus,  greater  overall  attenuation 
is  possible  if  n  is  made  an  even  number. 

From  experimental  tests,  the  attenuation  of  the  acous- 
tical  beam  per  reflection  is  on  the  order  of  10  dB,  at  a 
fused-silica  lead  interface.  Hence,  to  obtain  40  dB  atten¬ 
uation  4  reflections  are  required.  A  suitable  wedge  angle 
is  then  36  deg. 

The  delay  medium,  itself,  is  plate  or  bar  shaped  de¬ 
pending  upon  the  size  of  the  optical  aperture.  The  faces 
normal  to  the  optical  axis  are  ground  flat  and  parallel  to 
a  degree  specified  by  the  application.*  The  fused-silica 
material  from  which  the  delay  line  is  made  is  of  Schlieren 
quality  so  that  inhomogeneities  and  residual  strains  are 
kept  to  a  minimum.  These,  if  present  to  a  large  degree, 
can  adversely  affect  the  diffraction  pattern. 

The  size  of  the  optical  aperture  is  governed  by  two 
factors,  the  resolution  requirements  of  the  system  and  ul¬ 
trasonic  diffraction.  The  former  controls  the  integration 
time  T  of  the  modulator  where, 


D  being  the  dimension  of  the  aperture  in  the  direction  of 
elastic  wave  propagation.  The  integration  time  can  be 
shown  to  be  inversely  proportional  to  the  peak-to- first- 
null  width  in  the  diffraction  pattern  and, hence,  to  the 
resolution  of  the  electro-optical  system.  12  ""  » 

*  At  this  laboratory  the  severest  specifications  have  been 
l^Ox  flatness  over  any  two  inch  square  area,  and  face  paral¬ 
lelism  to  less  than  one  minute  of  arc.  These  specifications 
have  proved  adequate  for  the  processor  application. 
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Ultrasonic  diffraction  causes  the  elastic  bean,  to  spread 
appreciably  beyond  the  Fresnel  region.  In  this  circumstance, 
the  resulting  optical  diffraction  pattern  cannot  be  explained 
by  the  plane  wave  analysis  presented  here.  If  the  minimum 
width  of  the  optical  aperture,  and  thus,  the  width  of  the 
transducer,  is  made  large  enough,  however,  the  Fresnel  re¬ 
gion  can  be  extended  beyond  the  aperture  length  and  the  un¬ 
desirable  condition  eliminated.  Therefore,  the  minimum 
width  of  the  optical  aperture,  in  the  simple  light  modulator 
of  Fig.  8,  is  determined  by  the  requirement  for  maintaining 
a  Fresnel  diffraction  field  throughout  the  length  of  the 
aperture. * 

A  more  complex  modulator  configuration,  applicable  to 
real-time  processing  of  signals  received  by  a  radar  using  a 
planar  array  antenna,  is  shown  in  Fig.  9-  This  modulator 
is  capable  of  simultaneous  time  and  spatial  multiplex  modes 
of  operation.12  The  number  of  transducers  in  the  transducer 
array  is  governed  by  the  size  of  the  antenna,  and  constra  led 
by  the  ultrasonic  diffraction  effects  and  practical  limita¬ 
tions  of  the  optical  system.  The  multiple  absorbing  wedges 
are  employed  to  obtain  a  maximum  optical  aperture  in  a  given 

size  plate. 

In  order  to  maintain  proper  electrical  phasing  in  the 
transducer  array  (an  important  requirement  in  the  radar  pro¬ 
cessor  application),  a  common  crystal  serves  all  the  trans¬ 
ducers.  This  can  impose  extremely  severe  conditions  on  the 
fabrication  of  the  array,  especially  at  very  high  frequen¬ 
cies  where  the  crystal  will  be  very  thin  (  <  3°vi) ,  but  quite 

long  ( B  3  in . ) . 


#  it  is  assumed  that  adequate  light  intensity  is  available 
for  all  aperture  widths  so  that  the  diffraction  pattern  is 
detectable  within  the  system  constraints. 
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FIG  9  MULTI-TRANSDUCER  SOLID  LIGHT  MODULATOR  CONFIGURATION 
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In  both  the  simple  and  the  more  complex  modulator  con¬ 
figurations,  the  minimum  thickness  of  the  delay  line  will 
be  a-function  of  the  desired  modulation  index,  the  available 
driver  voltage  and  the  optical  specifications  for  the  plate. 
Precision  grinding  requires  that  the  plate  have  strong  flex¬ 
ural  resistance.  Hence,  even  though  the  strength  of  the 
elasto-optic  interaction  may  indicate  that  a  rather  thin 
plate  would  be  satisfactory,  the  fabrication  requirements 
may  actually  dictate  he  allowable  minimum  thickness. 

B.  transducer  considerations 

Quartz  transducers,  driven  at  their  fundamental  fre¬ 
quency,  will  be  used  to  generate  the  elastic  waves  in  the 
solid  light  modulators  under  discussion.  By  choosing  the 
proper  crystal  cut,  either  compressional  or  shear  waves 
will  be  produced.  The  X-cut  crystal  is  most  suitable  for 
the  generation  of  compressional  waves,  while  a  Y-or  AC-cut 
may  be  used  to  generate  shear  waves.  In  the  latter  case, 
however,  the  AC-cut  has  been  found  impractical  because  the 
thinness  of  the  crystal  at  100  MHz  frequency  makes  it  ex¬ 
tremely  difficult  to  handle.  The  treatment  to  follow, 
therefore,  will  be  concerned  with  X* and  Y-cut  transducers 

only. 

1.  general  Considerations 

The  resonant  frequency  of  an  X-or  Y-cut  crystal 
1.  inversely  proportional  to  the  thickneea  of  the  cryetal, 
according  to  the  formula, 
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where  v^  Is  the  velocity  of  propagation  of  the  elastic 
wave  in  the  crystal.  Given  the  resonant  frequency  in  MHz, 
tQ,  the  required  thicknesses  for  each  type  crystal  may  be 
calculated  from  the  relations j 

t  .  10*  x_cut 

fo 

m  Al  JkQf  Y-CUt 

fo 

with  tc  expressed  in  microns  (u).  Thus,  for  f^  •  100  MHz, 
the  X-cut  crystal  will  have  a  thickners  of  28. 5u  and  the 
Y-cut,  19. 2U. 

The  crystals  are  attached  to  the  fused-silica  de¬ 
lay  line  by  a  metallic  bond  of  either  lead  or  indium.  This 
bond  can  significantly  influence  th<*  transducer  performance, 
as  will  be  shown.  Where  a  50  P«r  cent  mechanical  bandwidth 
is  required,  the  free  crystal  face  may  be  left  unloaded. 

For  bandwidths  on  the  order  of  100  per  cent,  however,  load¬ 
ing  is  required,  with  a  consequent  reduction  in  force  at 
the  transducer-delay  line  interface. 

The  electromechanical  equivalent  circuit  for  an 
air-backed  X-or  Y-cut  transducer,  including  the  bond,  which 
is  assumed  lossless,  is  shown  in  Fig.  10, The  crystal 
parameters  are  the  static  capacitance  CQ,  the  electrome¬ 
chanical  transformer  ratio  <J>,  and  the  mechanical  impe¬ 
dance  of  the  crystal  Zc>  The  bond  parameters  are  the  me¬ 
chanical  impedance  of  the  bond  material  Z^,  and  the  wave¬ 
length  thickness  of  the  bond,  6.  The  element  Z^  repre- 
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FIG  10  ELECTROMECHANICAL  EQUIVALENT  CIRCUIT  FOR  X  AND  Y-CUT  QUARTZ 

transducers  including  thin  adhesive  bond 
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sents  the?  mechanical  impedance  of  the  fused-silica  delay 
medium,  which  is  real,  and  the  parameter  x  is  frequency 
dependent  and  defined  as 


X  £  7T  J- 

o 

where  f  is  the  electrical  excitation  frequency  in  Hz. 

In  quantitative  terms,  the  crystal  parameters,  as 
functions  of  the  resonant  frequency  and  the  effective  trans 
ducer  area  A  in  m2,  ares20 


1.425  x  104  Th 

Pf 

X-cut 

2.11  X  104  Ta 

Pf 

Y-cut 

60.6  f^A 

coulombs/m 

X-cut 

89.2  f^A 

coulombs/m 

Y-cut 

15.1  X  10®A 

kg/sec 

X-cut 

10.1  x  10°A 

kg/sec 

Y-cut 

The  mechanical  impedance  of  the  bond  and  the  de¬ 
lay  medium  depends  upon  the  type  of  wave  being  propagated 
as  well  as  the  choice  of  material.  When  compressional  waves 
are  excited, 


Zq  -  13. 1  X  10°A  kg/sec 

and,  according  to  a  relatively  recent  source,21 

2^  -  24.6  x  10°A  kg/sec 
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for  a  thin  lead  bond,  and 

Zb  -  18.3  x  10eA  kg/ sec 

for  a  thin  indium  bond.  The  corresponding  shear  wave  case 
is  covered  by  the  expressions: 


8.29  X  10®A 

kg/ sec 

8.0  X  10°A 

kg/ sec 

( lead) 

5.2  X  10°A 

kg/ sec 

( indium) 

where  Z.  (indium)  differs  from  the  expression  given  in 
the  source  cited  because  it  appears  to  be  incorrect  there 
on  the  basis  of  data  received  by  CUERL. 18 

The  wavelength  thickness  of  the  bond,  6,  is  a 
dimensionless  parameter  defined  by 


Where  t^  is  the  actual  bond  thickness  and  is  the 

wavelength  of  the  excited  wave  in  the  bond. 

The  values  of  the  frequency  dependent  parameters 
are  given  in  Table  I  for  the  two  frequencies  of  interest, 
80  and  100  MHz.  Determination  of  the  required  effective 
are*  depends  upon  other  considerations  Which  will  be  dis¬ 
cussed  subsequently. 
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TABLE  X 

X-j-  Y-CUT  TRANSDUCER  PARAMETER  VALUES  AT  80  AND  100  MHz 


PARAMETER 

X-Cut  (Compressional) 

Y-Cut  (Shear) 

80  MHz 

100  MHz 

80  MHz 

100  MHz 

t  Microns 
c 

35.6 

28.5 

24.0 

19.2 

C  A 

°  A  M* 

114  x  104 

142  x  104 

168.8  x  104 

211  X  10 4 

*/  Coulombs 
/  h 

M 

4848 

6060 

7136 

8920 

tc  *  thickness  of  transducer  crystal 
Cq  ■  static  capacitance  of  crystal 
$  *  electromechanical  transformer  ratio 
A  *  effective  transducer  area 
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2.  Transducer  Performance  with  Ze fOrTh lCkn£fi,£ 

When  the  bond  thickness  is  zero,  6  -  0  and  the  bond 
elements  in  the  equivalent  circuit  of  Fig.  10  have  no  effect 
upon  the  transducer  performance.  This  is  the  so-called 
"ideal"  case  which  is  most  important  in  determining  the 
nominal  characteristics  of  the  transducer  in  the  prelimin¬ 
ary  design  of  the  solid  light  modulator. 

The  force  F  at  the  transducer-delay  medium  in¬ 
terface,  in  this  circumstance,  is  given  by 

I  F_  .  _ i - 7  W 

I  Fo  (1  -  cot2  |)(1  +  c2tan2x) 


where 


F  *  2<J>V  , 
o 


(75' 


the  generated  force  at  f  -  f„  (i.e.,  at  the  crystal  reso 
nant  frequency),  and 


From  Eq.  ( 7 41,  the  frequency  response  of  an  X- 
and  Y-cut  transducer  radiating  into  a  semi-  inf  mite  f  used- 
silica  delay  medium  can  be  determined.  Figure  11  shows 
these  responses  as  a  function  of  the  normalized  frequency 
f/f  .  Note  that,  in  both  cases,  50  per  cent  bandwidths  are 
easily  achieved  without  the  aid  of  backing.  And, for  all 
practical  purposes,  the  attainable  unbacked  transducer  band 
width  is  the  same  for  both  transducers. 
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FIG  1 1  IDEAL  FREQUENCY  RESPONSE  OF  SINGLE,  AIR-BACKED  TRANSDUCERS 
RADIATING  INTO  NON-REFLECTIVE  FUSED-SILICA  MEDIUM 
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The  input  impedance,  at  resonance,  is  an  import¬ 
ant  parameter  because  it  determines  the  load  characteris¬ 
tics  for  the  transducer  driver.  Referring  to  the  equivalent 
circuit,  it  is  seen  that  when  f  •*  fQ,  the  only  load  on 
the  electromechanical  transformer  is  Z  .  Hence,  the  elec- 
trical  load  presented  to  the  driver  is  a  simple  parallel 
RC  network  whose  parameter  values  are  (neglecting  external 
factors) : 


Zq 

R .  =  — -a—  ohms 

1  4<t>2 

Ci  =  Co  pf 

when  Z  is  in  kg/sec  and  <t>  is  in  coulombs/m. 

The  average  power  absorbed  by  th*.  transducer  is 

given  by 

F  2 

Pav  ’  2^  ' 

when  the  excitation  is  sinusoidal  and  generates  a  force 
amplitude  F  .  At  the  resonant  frequency, 

F  *  2<J>V 
m  m 


by  Eq.  (75)  >  and  thus. 


2<t>2V  2 
m 


av 
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Since  both  C*  and  Zq  are  functions  of  the  area 
A,  the  parametric  dependence  of  the  average  absorbed  power 
may  be  expressed  more  specifically  as, 

P  -  56lAf  2V  2  x  10" 8  X-cut 

av  o  m 

«  1917Af  2V  2  x  10"®  Y-cut 
o  m 


where  P  is  in  watts. 

9V 

The  working  forms  of  the  equations  for  determining 

P  ,  R.  and  C,  at  80  and  100  MHz  are  given  in  Table  II. 
av'  l  i 

3.  Transducer  Performance  with  Non-Zero  Thickness 
Bond 

When  the  bond  thickness  is  not  zero,  but  some 
fraction  of  a  wavelength,  the  normalized  force  at  the  in¬ 
terface  may  be  expressed  in  the  form, 


F 

l|"(l  -  cot^)cos  2x6 

+  (^)cotf  sin  2x6 j 

(l 

1 

+  b= 

(’^)cofci  cos  2x6  " 

( 1  -  cot2^)sin  2x6  ?, 

where 


is  the  normalized  bond  impedance. 

At  the  crystal  resonant  frequency,  f  »  fQ,  the 
parameter  x  »  rr,  and  Eq.  (76)  reduces  to 
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TABLE  II 

X-&  Y-C*JT  TRANSDUCER  ELECTRICAL  PARAMETER  FORMULAS 
AT  RESONANT  FREQUENCIES  OF  80  AND  100  MHz  - 
ZERO -THICKNESS  BOND 


PARAMETER 

X-Cut  (Compressional) 

Y-Cut  (Shear) 

80  MHz 

100  MHz 

80  MHz 

100  MHz 

ohms 

,9.133 

A 

0.08Q 

A 

0.041 

A 

0. m 

A 

ci  ■  Co  Pf 

* 

O 

H 

X 

< 

H 

H 

142. 5A  x  104 

♦ 

2 

X 

& 

VO 

H 

* 

0 

H 

X 

2 

H 

CM 

Pav  watt. 

3.59AVm= 

5.61AV  2 
m 

12.27AV  2 
m 

!9.17AVm2 

■  input  resistance  of  transducer 
-  input  capacitance  of  transducer 
Pav  ■  average  absorbed  power 
A  -  effective  transducer  area  in  m2 
Vm  -  applied  voltage  amplitude 
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f  - 


(cos22tt6  +  b2«in227T6)~1  . 


This  normalized  force  is  plotted  in  Pig.  12  as  a  function 
of  6,  for  X-and  Y-cut  crystals  bonded  with  lead  and  in¬ 
dium  to  the  fused-silica  delay  medium.  It  should  be  noted 
that  the  difference  in  the  shapes  of  the  curves  is  solely 
dependent  upon  the  normalized  bond  impedance  b. 

The  LEAD-SHEAR  curve  exhibits  the  least  change  in 
|  F  |  Pq  |  with  6.  This  is  indicative  of  the  fact  that  the 
mechanical  impedance  of  the  lead  closely  matches  that  of 
the  fused-silica  in  shear  operation.  The  value  of  b  in 
this  case  is  0. 965. 

As  the  impedance  mismatch  increases,  the  variation 
of  |  F  |  Fq  |  with  changes  in  6  increases.  When  a  lead 
bond  is  used  with  an  X-cut  crystal,  b  *  1.88,  and  rapid 
variation  occurs,  reaching  a  minimum  of  F  -  0. 5  F  at 
6  ■  0.25  and  0.75.  With  an  indium  bond,  however,  the 
variation  is  significantly  smaller  since  b  ■  1.4.  It  ap¬ 
pears,  then,  on  the  basis  of  these  curves,  that  it  is  ad¬ 
vantageous  to  bond  an  X-cut  crystal  with  indium  rather  than 
lead  in  order  to  decrease  the  bond  thickness  sensitivity  of 
the  transducer.  In  both  cases,  as  will  be  shown,  a  wave¬ 
length  bond  thickness  less  than  6  *  0.025  will  be  neces¬ 
sary  if  maximum  response  and  bandwidth  are  to  be  achieved. 

The  INDIUM- SHEAR  curve  shows  the  sharpest  varia¬ 
tion  in  response  with  changes  in  6,  but  in  the  opposite 
direction.  This  is  characteristic  of  bonds  whose  b  values 
are  less  than  one  (here  b  -  0.627).  Outside  of  the  dif- 
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ficulty  of  bonding  with  lead,  there  ie  no  apparent  reaeon 
for  choosing  inc  um  ae  the  bonding  material  with  Y-cut 
transducers,  since  the  lead  bond  so  obviously  displays  an 
almost  ideal  character.  If  indium  is  used,  however,  maximum 
performance  will  only  be  achieved  if  6  is  again  less  than 

0.025. 

Although  all  the  curves  in  Fig.  12  are  symmetri¬ 
cal  about  6-0.5,  it  cannot  be  assumed  that  the  frequency 
response  will  also  be  symmetrical  about  this  value.  Quite 
the  contrary  turns  out  to  bo  the  case,  as  shown  below,  so 
that  any  attempt  to  obtain  optimal  transducer  performance 
must  be  made  in  the  direction  of  minimum  thickness  bonds. 

The  frequency  response  of  the  transducers  under 
consideration  can  now  be  obtained  from  Eq.  (76) ,  by  letting 
x  vary  over  the  range  §  <  x  <  ^  tor  .elected  value,  of 
6.  The  response  curve,  are  plotted  In  Figs.  13  “ 

pig,  13a  shows  the  variations  in  response  ioi  an 
X-cut  crystal  bonded  to  the  fused-siUca  delay  medium  with 
lead  of  thicknesses  6  -  0.01,  0.05  and  0.1.  The  6  -  0 
curve  is  also  plotted  to  establish  a  reference.  The  high 
sensitivity  to  bond  thickness  is  clearly  evident  along  with 
the  rapid  shift  in  peak  response  toward  lower  frequencies 
as  6  increases.  Reduced  bandwidth  and  increased  asymmetry 
accompanies  the  shift. 

The  formation  of  double-  and  triple-peaked  response# 
first  beginning  to  appear  in  the  6-0.1  curve,  is  shown 
clearly  in  Fig.  13b»  where  the  6  «  0.25#  °*5#  °*75  and  1,0 
curves  are  plotted.  The  6-0.5  curve  is  essentially 
single-peaked  but  its  bandwidth  is  significantly  less  than 
the  curves  for  6  <  0.25.  Thus,  its  possible  use  in  the 
light  modulators  of  concern  to  this  project  is  eliminated. 
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FIG.  13a  MECHANICAL  RESPONSE  OF  X-CUT  TRANSDUCER  WITH  LEAD  BOND 

TO  FUSED '  SILICA  DELAY  MEDIUM 
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HG  13b  MECHANICAL  RESPONSE  OF  X-CUT  TRANSDUCER  WITH  LEAD  BOND 

TO  FUSED -SILICA  DELAY  MEDIUM 
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If  operation  of  the  X-cut  transducer  with  lead 
bond  requires  full  50  per  cent  bandwidth,  centered  very 
close  to  the  crystal  resonant  frequency,  it  is  necessary, 
according  to  Fig.  13a,  that  the  bond  thickness  be  much  less 

than  6  •  0.05. 

As  a  measure  of  the  actual  thicknesses  involved, 
calculations  of  bond  thicknesses  in  microns,  at  80-  and  100-MHz 
frequencies,  have  been  made  and  are  presented  in  Table  III 
for  the  values  of  6  appearing  in  the  response  curves.  It 
can  now  be  noted  that  at  100  MHz,  is  only  1.  In  when 

6  -0.05  with  lead. 

The  mechanical  response  of  a  Y-cut  crystal  with 
lead  bond  is  shown  in  Figs.  14a  and  b.  The  good  mechanical 
match  between  the  bond  and  the  medium  maintains  almost  sym¬ 
metrical  response  about  f  -  fQ»  regardless  of  the  value  of 
6.  Bandwidth  is  also  preserved  in  all  cases.  Hence,  the 
bond  thickness  has  negligible  effect  upon  the  performance 
of  this  transducer.  This  offers  a  distinct  advantage  when 
spatially  multiplexed  modulators,  Which  require  large  com¬ 
mon  crystals  uniformly  bonded  to  the  delay  medium,  are  under 

consideration. 

Improved  congressional  wave  transducer  response 
is  possible  when  indium  is  substituted  for  lead  in  the  bond. 
This  is  clearly  indicated  in  Figs.  15«  »nd  b,  and  reinforces 
the  prediction  of  this  result  made  on  the  basis  of  the  curves 
in  Fig.  12.  Indium  is  easier  to  bond  with  and,  for  a  given 
6,  has  a  slightly  greater  thickness  (Table  III). 

Figures  l6a  and  b  show  the  response  to  be  expected 
from  an  indium-bonded  Y-cut  crystal.  Because  b  <  1,  the 
curves  for  small  6  tend  to  have  peak  responses  at  frequen- 
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FIG  14b  MECHANICAL  RESPONSE  OF  Y-CUT  TRANSDUCER  WITH  LEAD  B  >NL 

TO  FUSED  'SILICA  DELAY  MEDIUM 
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FIG  15b  MECHANICAL  RESPONSE  OF  X-CUT  TRANSDUCER  WITH  INDIUM  BOND 

TO  FUSED -SILICA  DELAY  MEDIUM 
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FIG.  160  MECHANICAL  RESPONSE  OF  Y-CUT  TRANSDUCER  WITH  INDIUM  BOND 

TO  FUSED  -  SILICA  DELAY  MEDIUM 
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F,G  16b  MECHANICAL  RESPONSE  OF  Y-CUT  TRANSDUCER  WITH  INDIUM  BOND 

TO  FUSED -SILICA  DELAY  MEDIUM 
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cies  above  fQ.  However,  when  6-0.25  and  above,  the 
response  becomes  symmetrical  about  f  -  fQ#  with  generally 
reduced  bandwidth.  It  is  of  interest  to  note  that  at 
6  -  0.25,  F  is  approximately  60  per  cent  greater  than  it 
would  be  in  the  ideal  case,  6-0.  This  enhancement  in 
the  response  might  be  useful  in  certain  narrow-band  appli¬ 
cations,  but  is  of  little  concern  here.  Table  III  indicates 
that  at  100  MHz,  -  1.775m-  when  6  -  0.25.  Therefore, 
it  is  easy  to  deposit  an  indium  layer  of  this  thickness 
with  present  deposition  techniques. 

c.  MODULATION  INDEX  VS  APPLIED  VOI/TAGE 

The  relationship  between  the  modulation  index  and  the 
voltage  applied  to  the  transducer  can  now  be  established 
by  Eqs.  (63),  (64)  and  (75).  It  is  assumed  that  the  exci¬ 
tation  is  sinusoidal,  with  amplitude  Vm  and  frequency  fQ, 
and  the  transducer  bond  approximates  the  ideal  zero  thick¬ 
ness.  Then,  by  Eq.  (75)# 


m 


2<W. 


m 


and  the  maximum  stress  Tm  is  just 


2(T)V, 


m 


(77) 


where  it  is  further  assumed  that  the  transducer  action  is 
uniform  over  the  effective  area  A. 

Applying  Eq.  (77)  to  Eqs.  (63)  and  (64)  gives  the  fol 
lowing  relationships  for  L  and  T  waves,  respectively! 
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*  «  2qL(^)V 

*m  _  \  E  '  o  H  v  A  '  m 

lL»  L 


*  -  2pL(““)v 

hn  '  E  '  o  *  '  A  '  m 

2L  L 


m  _  in 
lT  i 


ff  “<£>*. 

L 


(78) 


(79) 


where  $  is  the  electromechanical  transformer  ratio  for 
x 

the  X-cut  transducer  and  <J>  is  that  for  the  Y-cut  trans- 

y 

ducer.  All  other  terms  correspond  to  those  defined  pre¬ 
viously,  and, since  only  magnitudes  are  of  interest  here, 

the  signs  of  the  #  have  been  neglected. 

m 

The  working  forms  of  these  equations,  at  80-and  100-MHz 

resonant  frequencies,  can  be  obtained  by  substitution  of 

the  appropriate  values  for  all  the  quantities  on  the  right 

except  L  and  V  ,  which  are  independent  parameters.  The 

m 

appropriate  values  are  listed  in  Tables  I  and  IV,  where 
the  value  for  is  taken  as  that  of  the  He-Ne  laser 

ll 

light  source  being  used  in  the  experimental  program  at 
CUERL.  with  these  values,  Eqs.  (78)  and  (79)  reduce  tos 

f  =80  MHz 
o _ 

=  27.6  X  10"*L  V  (80) 

m  m  ' 

lb 
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TABLE  IV 

VALUES  OF  CONSTANT  PARAMETERS 
USED  IN  DERIVATION  OF  EOS.  (80)  &  (Blj 


\  **  0.633  X  10"®  m 

L 

0  c:  0.  164 

E  =  0.76  x  10 11  newtons/m2 


n  =1. 46 
o 


q  *  0.088 

p  *=  0. 197 


Measured  with  0.589u 
light  wavelength* 


B  =  3.51 


x  10"l2m2/newton 

measured  with  0.644u. 
light  wavelength* 


♦Primak,  W. ,  and  Post,  D. ,  "Photoeia.tic  Constants  of 
Vitreous  Silica  and  Its  Elastic  Coefficient  of  Re 
tractive  Index,"  Journal  of  Applied  Physics,  Vol.  3  » 

No.  5»  May  1959. 


T- 3/321 


-78- 


COLUMBIA  UNIVERSITY— ELECTRONICS  RESEARCH  LABORATORIES 


f  =  80  MHz  (Cont'd. ) 


=  61.7  x  10“4L  V 

Xu 

^  ■  25  x  io~4l  v 

in  _  in 

2? 


f  B  100  MHz 

o  _  _ 


=  3^.5  x  10“4l  v 

m  xn 

aL 

if  =  77.2  x  io”4l  v 

m  xn 

2L 

iff  b  ip  e  31.2  X  10  4L  V 
m  -  m  _ 
lT  2T 


(80) 


(81) 


where  L  is  in  cm  for  convenience,  and  V  is  in  volts. 

in 

As  am  example  of  the  voltages  required  to  achieve  a 

desired  modulation  index  f  <0.3  radian,*  the  second  and 

m  _ 

third  formulas  in  Eqs.  (80)  and  (8l)  are  plotted  in  Figs.  17 
and  18,  for  selected  values  of  L,  where  the  direction  of 
polarization  of  the  incident  light  has  been  chosen  to  be  in 


*  Lambert12  has  shown  that  to  maintain  a  dynamic  range  of 
>  40  dB,  muat  be  <0.3  radian.  The  dynamic  range  is 

defined  as  the  ratio  of  the  peak  signal  intensity  to  the 
peak  spurious  response  intensity  in  the  image  plane  of  the 
electro-optical  processor. 
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FIG. 1 7  MODULATION  INDEX  VS.  COMPRESSIONAL  TRANSDUCER  VOLTAGE 
WITH  L  AND  f  AS  INDEPENDENT  PARAMETERS 


MAXIMUM  VOLTAGE  ACROSS  TRANSDUCER  (VOLTS) 


FIG  18  MODULATION  INDEX  v*.  SHEAR  TRANSDUCER  VOLTAGE  WITH 
L  AND  f  AS  INDEPENDENT  PARAMETERS  . 
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the  2 -direction,  i.e.,  normal  to  the  direction  of  propaga¬ 
tion  of  the  elastic  waves. 

It  is  seen  that  for  a  given  modulation  index,  the  re¬ 
quired  voltage  in  the  L  wave  mode  of  operation  is  signif¬ 
icantly  less  than  in  the  T  wave  mode.  This  can  be  ex¬ 
plained,  theoretically,  by  considering  Eqs.  (61),  (78), 

and  (79). 

From  Eq.  (61),  the  stress-optic  constant  is  related 
to  the  Neumann  strain-optic  constants  through  the  equation. 


B 


‘(q  -  P)  -  '  1)P  * 


Substituting  the  numerical  values  for  p 


and  q  then  crives 


b  I  -  0.55M 


l 


)no*P 


This  quantity  can  now  be  compared  to  the  corresponding  quan 
tity  in  the  formula  for  *m  .  If 


bl 


1+0 


)"oEp 


then 


!  B 


0.277  Bl  . 
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i 

Hence,  even  though  <t>  <  <J>  ,  for  sny  given  values  of  fQ 

and  A,  the  above  relationship  indicates  that  the  L  wave 
interaction  will  be  markedly  stronger  than  the  T  wave 
interaction . 

The  high  voltages  required  to  achieve  the  limiting 

modulation  index,  f  *=0.3  radian,  in  Fig.  18  shows  that 

m 

the  shear  mode  of  operation  is  highly  inefficient  and,  ac¬ 
cording  to  the  power  formulas  ir.  Table  II,  can  be  expected 
to  absorb  a  relatively  large  amount  of  pi  ’er,  even  when  the 
transducer  area  is  small. 

Considerable  improvement  in  this  condition  can  be  ob¬ 
tained  by  operating  ir  the  compressional  mode  where,  by 
reference  to  Fig.  17,  the  required  voltage  for  a  given  ^ 
an <3  l  is  seen  to  be  less  than  one-half  that  for  the  com¬ 
parable  shear  mode. 

*  i 

D.  INTERNAL  REFRACTION  EFFECTS 

So  far,  in  this  study  of  a  solid  light  modulator,  the 
effect  of  a  refractive  index  gradient  in  the  elastic  wave 
field  has  been  neglected.  As  a  result,  it  has  been  shown 
that  the  light  is  only  phase-modi lated  after  passing  through 
the  field  and  the  diffraction  can  be  explained  by  the  simple 
Raman -Nath  theory. 

Actually,  the  existence  of  refractive  index  gradients 
in  the  sinusoidally  varying  elastic  field  will  tend,  through 
/  refraction,  to  bunch  and  spread  apart  the  light  rays  with  a 

f 

periodicity  equal  to  that  of  the  field.  Since  the  emerging 
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light  intensity  is  directly  related  to  the  ray  density,  am¬ 
plitude  modulation  will  occur  simultaneously  with  the  phase 
modulation.  If  this  amplitude  modulation  is  large  enough, 
it  can  adversely  effect  the  simple  relationship  between  the 
diffracted  light  amplitude  and  the  modulation  index.  There¬ 
fore,  identification  of  the  controlling  parameters  is  im¬ 
portant  in  order  to  establish  the  necessary  criteria  for 
limiting  the  modulation  to  an  acceptable  level. 

Analytical  investigations  of  internal  refraction  ef¬ 
fects  have  proceeded  along  two  lines  of  approach.  In  one, 
the  geometrical  problem  of  light  rays  bending  in  a  variable 
refractive  index  field  has  been  treated.  In  the  other,  the 
effect  of  such  bending  on  the  diffraction  pattern  itself 
was  of  primary  concern. 

The  former  approach  yields  a  picture  of  the  bending 
over  an  elastic  wavelength  as  a  function  of  known  parameters. 
Therefore,  the  extent  of  the  bunching,  for  a  specific  set 
of  parameter  values,  is  made  clear  and  a  limiting  criteria 
for  simple  Raman-Nath  diffraction  can  be  established. 

The  latter  approach  shows  that,  for  normal  light  inci¬ 
dence,  the  intensity  of  the  diffraction  fringes,  but  not 
their  form,  will  be  affected  as  a  result  of  the  refraction 
phenomenon.  A  second,  stronger,  criterion  for  limiting  the 
amplitude-modulation  to  an  acceptable  level  can  then  be  es¬ 
tablished  by  noting  the  variation  in  the  first-order  inten¬ 
sity  from  that  predicted  by  the  simple  Raman-Nath  theory, 
as  the  governing  parameters  are  changed. 
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The  geometrical  theory  of  internal  refraction  derives 
trom  the  work  of  Lucas  and  Bicquard. ‘ *  Although  their  treat- 
ment  is  quite  general,  a  solution  for  the  trajectory  equa¬ 
tion  was  obtained  only  for  the  case  where: 


1.  The  index  of  refraction  could  be  represented  by 
the  equation. 


27T 

n  =  n  +  6n  cos  - —  x  . 
o  m  X 

e 


2.  6n  ((  n 
m  o 


Since  these  conditions  apply  in  the  solid  light  modulator 
under  consideration,  the  results  obtained  from  their  theory 
should  be  valid. 

The  computed  light  ray  trajectories  over  an  elastic 
wavelength,  as  a  function  of  the  parameter23 


M 


_1_ 

X 


V 


2tt^  X.L 

m  l 


n 


(82) 


is  shown  in  Fig.  19.  The  first  crossing  point  occurs  at 

7T 

M  =  1,57  c  —  .  The  significance  of  this  crossing  point  is 
that,  if  the  parametric  values  in  Eq.  (82)  are  such  that  a 
value  of  M  ,>  1.57  is  obtained,  the  light  emerging  from 
the  modulation  will  appear  to  come  from  a  series  of  line 
sources  spaced  Xfi  apart.  Thus,  the  nature  of  the  diffrac- 
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FIG.  19  LIGHT  RAY  TRAJECTORIES 
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tion  will  change  significantly  as  the  modulator  is  trans¬ 


formed  from  a  pure  phase  grating  into  a  line  and  phase 


grating. 

If  M  =  1.57  is  now  considered  a  limiting  criterion 
for  simple  Raman-Nath  diffraction,  the  critical  relation¬ 
ship  between  the  modulation  index  and  L  can  be  determined 
from  Eq.  (82),  with  M  -  -  .  Hence, 


n  *  2 
o  e 


(83) 


This  expression  is  Willard's  criterion24  for  simple  Raman- 
Nath  diffraction,  in  slightly  modified  form.  It  indicates 
the  maximum  modulation  index  that  will  still  maintain  the 

simple  Raman-Nath  diffraction  for  a  given  L,  once  the 

» 

other  parameter  values  are  known. 

Equation  (83)  is  plotted  in  Figs.  20  and  21  for  L 
and  T  wave  modes  of  operation,  respectively.  In  each 
plot,  curves  for  the  80  and  100  MHz  frequency  are  given. 

The  modification  to  the  first-order  peak  intensity  of 
the  diffraction  pattern  due  to  internal  refraction  has  been 
derived,  in  a  simple  way,  by  Rao  and  Murty.25  For  normally 
incident  l^ht  and  <[  1  their  expression  can  be  written 
in  the  form. 
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FIG  20  LIMITING  VALUES  OF  MODULATION  INDEX  FOR  NORMAL  DIFFRACTION 
VS  WIDTH  OF  ELASTIC  WAVE  FIELD  L,  FOR  LONGITUDINAL  WAVES 
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BASED  ON  WILLARD  S  CRITERION 
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FIG.  21  LIMITING  VALUES  OF  MODULATION  INDEX  FOR  NORMAL  DIFFRACTION 
VS.  WIDTH  OF  ELASTIC  WAVE  FIELD  L,  FOR  SHEAR  WAVES 
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(84) 


(85) 


where  Iim(LAe2)  i8  the  niaxiraum  first-order  intensity 

when  internal  refraction  is  taken  into  account,  and  I  (0) 

in* 

is  the  intensity  given  by  the  simple  Rare an -Nath  theory, 
e, , 


I 


(86) 


when  <  1. 
rm  N 

Figure  22  is  a  plot  of  I  (L/X  2)/l  (0)  vs  L  for 

2,m  e  jin 

80  and  100  MHz  in  the  shear  mode  of  operation.  Note  that 
severe  reductions  in  the  first-order  intensity  will  occur 
for  the  larger  transducers  even  though  Willard's  criterion 
is  satisfied. 

The  validity  of  this  plot  can  be  tested  by  applying 
tho  most  severe  practical  conditions  to  the  inequality, 

Eq.  (85).  Letting  X  =  4.7  x  10_3cm  (80  MHz),  L  =  0.2  cm, 

G 

and  Tp  -  0. 3, 
ro 
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FIG.  22  FIRST-ORDER  INTENSITY  AFTER  INTERNAL  REFRACTION 

BY  SHEAR  WAVES 
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Wv*"  °-l22«  1  • 

L 


Hence,  Eq.  (84)  is  valid  over  the  entire  useful  range  of  L. 


Hie  corresponding  plot  for  the  longitudinal  wave  mode 
of  operation  is  shown  in  Fig.  23.  It  is  immediately  evident 
that  a  significantly  smaller  reduction  in  intensity  will 
take  place,  for  a  given  L,  with  this  mode  of  operation 
when  compared  to  the  shear  mode.  It  should  be  noted,  how¬ 
ever,  that  \he  condition,  Eq.  (85)#  is  violated  when 
L  <  0.2  cm  at  80  MHz.  This  presents  no  problem  because, 
from  a  design  viewpoint,  the  L  wave  transducers  will  al¬ 
ways  call  for  larger  values  of  L  as  a  result  of  the  opti¬ 
mal  design  considerations  to  be  discussed  in  the  next  sec¬ 


tion. 


The  expression  for  the  first-order  peak  intensity  can 
be  written  in  terms  of  the  peak  intensity  of  the  zero-order 
by  using  Eq.  (86).  Hence, 


"V 


2n  X  2 

o  e 


nXLL 


(87) 


2n  X  2 
o  e 


where  it  is  seen  that  the  effect  of  internal  refraction  is 
to  "weight”  the  intensity  obtained  from  the  simpler  theory 
without  changing  the  form  of  the  diffraction  pattern. 


Ai 
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FIG  23  FIRST-ORDER  INTENSITY  AFTER  INTERNAL  REFRACTION 
BY  LONGITUDINAL  WAVES 
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E.  qptti«t7ATTCM  OF  THE  TRANSDUCER  WIDTH 

the  transducer  width  is  equal  to  the  width  of  th 
elastic  wave  field,  L,  when  beam  spreading  is  constrained 
by  limiting  the  length  of  the  elastic  wave  propagation  in 
the  modulator  to  the  Fresnel  region. 

The  effect  of  this  parameter  on  the  first-order  peak 
intensity  ha.  been  indicated  in  the  previous  discussion, 
specifically  by  Eqs.  (62)  and  (87),  where  it  is  seen  that 
increasing  L  will  increase  but,  at  the  same  time, 

decrease  the  value  of  the  weighting  function  at  a  faster 

rate  Obviously,  some  optimal  value  of  L  must  exist 

.i„  in  Ea  (87)  is  maximized.  To  ob- 
which  the  intensity  ratio  m  Eq.  V  () 

tain  this  optimal  value,  consider  the  effective  modulation 
index  defined  as 


—  sm  v  . 

ip  =  -  t 

v 


m 


where 


7TXL 

o  e 


with  respect  to  L  and  setting 
the  derivative  equal  to  zero  gives. 


Differentiating  i>m 


dip 


m 


dL 


[cos  v  sin  v  1  dy.  ^ 
~  v2  J  dL  Vm 


,  d^ 
sin  v  _ 15.  _  n 

+  dL  - 
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or,  upon  substituting  for  v  and  > 


— —  =  cos  v  =  0  . 


Hence, 


v 


7r  2HL 

2*2*' 


Since  v  >  7r  would  carry  tin  v/v  out  of  the  main 
lobe,  the  only  satisfactory  solution  is  v  *=  ”  Thus, 


n  *  * 

o  e 


(88) 


This  expression  must  now  be  tested  to  see  that  both 
Willard^  criterion,  Eq.  ^83)*  and  the  condition  given  in 
Eq.  (85)  will  not  be  violated  for  modulation  indices,  y^, 
<0.3  radian.  Substituting  Eq.  (88)  into  Eq.  (83),  we  have 


^  noV 

8XLLopt 


=  0.39  radian 


which  implies  no  restriction  in  this  regard.  Substituting 
Eq.  (88)  into  Eq.  (85)  then  gives. 
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_m 

2n 


«  1 


7T 

as  a  second  criterion  on  ip  •  With  ip  =  x  ,  the  left 
^  m  mo 

side  equals  or  0.0625  «  1.  Hence,  Lopt  a 

valid  design  parameter  which  will  assure  maximum  electro¬ 
optic  performance  in  the  modulator. 

No  consideration  has  been  given  to  elastic  wave  at¬ 
tenuation  (which  would  supply  a  multiplicative  factor  to 
ip  )  because  at  the  frequencies  and  optical  apertures  of 
interest  the  effect  of  this  attenuation  in  fused-ti lica 
is  negligible. 20  Also,  since  the  attenuation  factor  is  in¬ 
dependent  of  L,  it  would  have  no  effect  upon  the  expres¬ 
sion  for  L  .  given  in  £q,  (88). 

opt 

The  weighting  function  in  Eqs.  (84)  and  (87)  can  now 

be  evaluated  and  L  ^  determined  for  L  and  T  wave 

opt 

modes  of  operation  at  80  and  100  MHz.  In  all  cases  it  is 

2 

seen  that  sin  v/v  *=  —  =  0.  6j7.  Therefore, 

7T 

w.p.  -  j4^]2  -  0.406. 

Referring  to  Figs.  22  and  23.  the  values  of  LQpt  are 
then : 

L  =  1.3  cm  (80  MHz,  L  mode) 
opt 

*  0.83  cm  (100  MHz,  L  mode' 

=  0.5  cm  (80  MHz,  T  mode) 

=  O.33  cm  (100  MHz,  T  mode). 
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The  theory  of  birefringence  in  a  transparent,  amorphous  solid  is 
developed  from  fundamental  principles  in  order  to  obtain  a  design  for¬ 
mulation  for  spatial  light  modulators  that  will  operate  at  100  MHz 
center  frequency  with  at  least  50  per  cent  bandwidth.  Both  longitu¬ 
dinal  and  shear  wave  elasto-op* ical  interactions  are  treated. 

Relations  for  the  phase  modulation  index  and  the  diffracted  light 
intensity  and  polarization,  as  functions  of  the  modulator  parameters 
®r*  presented  for  the  case  of  normal  light  incidence.  The  effect  of* 
internal  refraction  on  the  first-order  peak  intensity  is  then  quanti¬ 
tatively  examined.  From  this  consideration,  a  criterion  for  obtaining 
maximum  optical  performance  in  the  modulator  is  derived. 

An  analysis  of  the  quartz  transducers  used  to  generate  the  elas¬ 
tic  waves  in  the  modulators  is  presented  as  an  integral  part  of  the 

design  formulation.  The  effect  of  lead  and  indium  bonds  on  transducer 
response  is  treated  in  detail. 
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